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Q I Abstract 
CJ ■ 

By means of rather general arguments, based on an approach due to Derrida 
that makes use of samples of finite size, we analyse the effective diffusivity and 
^ I drift tensors in certain types of random medium in which the motion of the 

■ particles is controlled by molecular diffusion and a local flow field with known 

statistical properties. 

The power of the Derrida method is that it uses the equilibrium probability 
distribution, that exists for each finite sample, to compute asymptotic behaviour 
at large times in the infinite medium. In certain cases, where this equilibrium 
situation is associated with a vanishing microcurrent, our results demonstrate 
the equality of the renormalization processes for the effective drift and diffusivity 
tensors. This establishes, for those Ward identity previously verified only 

to two-loop order in perturbation theory in certain models. 

The technique can be applied also to media in which the diffusivity exhibits 
spatial fluctuations. We derive a simple relationship between the effective diffu- 
sivity in this case and that for an associated gradient drift problem that provides 
an interesting constraint on previously conjectured results. 
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1 Introduction 



Problems of flow and diffusion in random environments are of great interest for many 
reasons both mathematical and physical Jl], §]. Possible applications range from 
turbulent diffusion to ionic and molecular diffusion in the presence of trapped 



ions, dislocations and other impurities [|T0 



A typical continuum model, that represents some aspects of these physical systems, 
comprises a particle in a D dimensional space subject to molecular diffusion, k^, to- 
gether with a local drift velocity, u(x) . The motion of a particle at position x is given 
by 

i = u(x)+w(t) , (1) 
where w(t) is a white noise process that satisfies 

{wi{t))=0 and {wi{t)wj{t')) = 2Kij6{t - 1') . (2) 

The flow field u(x) is time independent and assumed to have certain statistical proper- 
ties, the details of which depend on the particular model. However we do assume that 
the system is homogeneous so that, for example, 

(Mi(x)nj-(x')) = Aij(x - x') , (3) 

and (u(x)) is independent of position. Frequently the statistics of the velocity field are 
assumed (for simplicity) to be Gaussian in character so that only the mean velocity and 
its two-point correlator are required to determine all correlators. We use angle brackets 
to indicate an average over the white noise ensemble and/or over the ensemble of 
samples of the medium as appropriate. The intended version of the averaging procedure 
will be clear from the context. 

The computational problem we address then, is to use the above statistical information 
on the model to deduce the effective drift and diffusivity that control the motion of the 
particles at large distances and times - the bulk properties of the medium. There has 



been a great deal of work on this problem in|-|T6[. In this paper we develop a 

continuum version of an approach by Derrida [l^ that he applied to lattice models 
in order to elucidate relationships that may exist between the effective long range 
parameters of the theory, one example of which is the "Einstein relation" between drift 
and diffusivity p|, ^. 
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2 Derrida Argument for the Continuum 



The idea behind Derrida's argument []T^, which is drawn from statistical physics, is 



that the bulk properties of an infinite medium can be captured in a sample of finite size 
provided this size is much larger than the typical correlation lengths of the statistical 
fluctuations of the velocity field u(x) . The size of the sample provides an infra-red 
cut-off, to use the language of field theory. The assumption being made is that the 
effective parameters are insensitive to the infra-red cut-off in a limit in which it becomes 
very large. Of course the presence of important long range correlations would lead to a 
dependence of the effective parameters on the sample size (and geometry). These are 
circumstances in which we would expect to see anomalous diffusion 
Another important issue in relation to the removal of the infra-red cut-off is the extent 
to which this limit commutes with the large time limit invoked in the definition of 
bulk quantities |l|] . We do not address this point mathematically in this paper but 
rely, for the moment, on physical intuition. The commutation of these limits is a very 
important issue both from a mathematical and a physical point of view and should be 
examined more thoroughly. 

We will assume, therefore, that the medium is represented by a velocity field Ui(x) 
that is defined over a basis region Q and that space is tessellated by this region under 
displacements I that lie in a Bravais lattice B . For the purposes of visualization it is 
sufficient to suppose that Q is cube of side L . The velocity field Uj(x) then satisfies 



and we assume that fl which contains the origin, is very large relative to the basic 
correlation length in Ui(x). 

A blob of particles is released from the origin at t = . They move independently 
according to eq(|l|) and are free to cross the boundaries of fl and it's replicates. The 
resulting probability distribution, P(x, t), satisfies 



All the moments of the distribution can be obtained from the moment generating 
function of the distribution 



Mj(x + 1) = Mi(x) , 



d_ 

di 



P(X, t) = di{Kijdj - Mi(x))P(x,t) 



(4) 




(5) 
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The technical step suggested by Derrida is to re-write this integral over all space as a 
sum over replicates of the basis region Q, thus 

Z(k)= / d^xye-*-("+')p(x + l,t) . (6) 
Jn J 

Z(k)= [ d^xW{ic,t)= f d^xe-'^-''R{x,t) , (7) 

i?(k,x,t) = ^e-^''-'P(x + l,t) , (8) 



That is 



where 



and 

W(k,K,t) ^ e-'^-''R{k,yi,t) , (9) 
The probability distribution can be reconstructed as 

(0e--Z(k) . (10) 

Note that there is no restriction on the integration range of k . 

Making use of the periodicity properties of the velocity field we can show that 

d 

— i?(k, X, t) = di{Kijdj - Ui{x))R{k, x, t) , (11) 
where i?(x, t) satisfies the boundary condition 

R{k, X + a, t) = e'^ ^'Rik, x, t) , (12) 
with a an element of the Bravais lattice. 

As i — > oo we expect that time dependence will disappear. In that limit 

di{Kijdj - Ui{x.))R{k, X, t) = . (13) 
This imphes that for large t 

i?(k,x,f) oc Qo(x) (14) 

where Qo^x) is the static probability distribution in Q . The advantage of the Derrida 
method is precisely that one may make use of this static equilibrium distribution to 
evaluate asymptotic behaviour. 

However Qo{yi) is periodic under displacements in the Bravais lattice and, except for 
the case k = 0, R{k, x, t) is not periodic. Hence we must have 

P(0,x,i)^Qo(x) , (15) 



and 

R{k,x,t) . (16) 

when k 7^ . All of this is entirely consistent with continuity in k at finite time. 
The equation satisfied by W(k,yi,t), which is periodic on the Bravais lattice, is 

|»'(k,x,i) = « + *)(K«(a, + .*,)-«.(x))iy(k,x.«) . (17) 

When k = we know that the operator on the right has an eigenfunction with a zero 
eigenvalue, namely Qo(x), the stationary probability distribution. When k 7^ this 
eigenfunction is perturbed and the eigenvalue is moved away from zero. We denote the 
eigenfunction by Q{x.) and the eigenvalue by — /i(k) so that yu(0) = and 

(9, + iki){K,,{dj + ikj) - M.(x))Q(x) = -/i(k)Q(x) . (18) 

An examination of the (trivial) case with constant velocity field suggests that near 
k = there is a gap in the spectrum of this operator of magnitude ~ (27r/L)^ provided 
|k| much less than 2tt/L . We assume this to be true and conclude that up to expo- 
nentially damped corrections that are 0(e~''^'^/^^^*), the solution of eq(0) in which we 
are interested is 

W{k,x,t) = g(x)e-'^('^)* . (19) 
We now impose the normalization condition 

/ c/^xg(x) = 1 , (20) 

to obtain the resulting formula for the moment generating function is 

Z{k) = e-f"^^^' , (21) 

where we have taken into account eq(^) . The corresponding formula for the proba- 
bility distribution in the unbounded (but tessellated) space is therefore 

^(x,t) = / (0je^'''""^'^* • (22) 

For large t the asymptotic form of P(x, t) is determined by the terms in /^(k) up to 
O(k^) . The remaining terms of higher order determine the approach to this asymptotic 
form in inverse powers of ^/t in the general case. 

It should be noted that on the basis of the above argument the asymptotic behaviour 
is guaranteed to set in only when t is so large that only values of |k| much less than 
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271 / L yield important contributions to the integral in eq(p2D . This means that we are 
computing asymptotic behaviour that becomes apparent after the probability distri- 
bution of the particle has spread over more than one replicate of . It is the physical 
assumption of the Derrida method that, when the samples are sufficiently large and 
there is no anomalous diffusive behaviour, this asymptotic behaviour sets in effectively 
much earlier and therefore there is no change to the estimated bulk diffusivity and drift 
when the infra-red cutoff is removed. Proving this mathematically remains to be done. 
Any such proof should also tell us about the physically important issue of how effects 
on different scales in the problem interact with one another. 



3 Asymptotic Behaviour 

We can calculate these terms in /i(k) by perturbation theory. Set 

g(x) = go(x) + gi(x) + g2(x)-f--- , (23) 

and 

/i(k) =^i(k)+/i2(k) + --- , (24) 

where the suffix on a term indicates its order in powers of k . Note that there is no 
term /iq and that the normalization of go(x) means that 

j rf^xgi(x) = j d^xg2(x) = . (25) 

This means that Qi, Q2 and all higher corrections lie in what we will refer to as the zero 
weight subspace. We can obtain the equations for the higher corrections by expanding 
both sides of eq(p!8D in powers of k . We find 

dl{Kimdm - M;(x))gi(x) + iki{2Kimdm " Mi(x))go(x) = -/ii(k)go(x) , (26) 

and 

dl {nimdm - (x) ) g2 (x) +iki {2Kimdm -Ul{y.))Qi (x) - ki km HlmQo (x) 

(27) 

= -/i2(k)go(x) - /ii(k)gi(x) . 

For completeness we note that the general recursion formula {n ^ 2) is 

di{Kimd.m-Ui (x))g„(x) +iki {2Kimdm-Ui (x))g„_i (x) (28) 
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= -/i„(k)Qo(x) (/W2(k) - kikmHim)Qn-2{^) - /il (k)Q„_i (x) . 

After integration over the base cell fl, eg (p6D yields 

yLii(k) = iki / (i^xM;(x)Qo(x) = ifcjMz (29) 

where -u^ is the average value of the velocity in the stationary probability distribution 
Qo{yi) . Note that the vanishing, under integration, of the gradient terms is assured by 
the periodicity and continuity of the random field. 

To compute the next correction we need the Green's function (^(x, x') that is the inverse 
of the operator 

H = di{Kimdm - Ui{x)) 

in the zero weight subspace. The existence of this Green's function depends on our 
physically plausible assumption that the kernel of H is one-dimensional and is therefore 
spanned by (5o(x) . 
We have then 

diinir^d^m - «Kx))G(x,x') = -5(x,x') = -5(x - x') + , (30) 

where V is the volume of 1) . It follows then that 

gi(x) = tki f t/^x'G(x, x'){2t,imd'^ - (m^x ) - ui))Qo{x') (31) 
Jn 

The evaluation of fi2 is achieved by integrating eq(27) with respect to x over Q . We 
have 

1^-2 = kiKimkm + iki / d^xui{x.)Qi{x.) (32) 
Jn 

That is 

1^2 = kiKirakm - hkm / d-^xd-^x'^^ (x)G'(x, x') (2k„„9^ - (um(x') - Um))Qo(x') (33) 

Jn 

These results tell us that the mean velocity of the cloud is Ui and the effective diffusivity 
tensor is 



i ||^d^xd^x'M,(x)G(x,x')(2/€„„9; - M„(x'))go(x') + (/ ^ m)| (34) 
where we have set 

ui{x) = ui{x) - ui{x) , (35) 

and used the result 

/ d^xG(x, x') = . (36) 
Jn 
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4 Gradient Flow 



A case of particular interest is that of gradient flow where the local drift depends 
linearly on the gradient of a scalar field. Thus 

Ui{x) = A;„9^0(x) . (37) 

In this case the equilibrium microcurrent is 

M^) = -{l^lmdm - \lradm<P{y^))Qo{y^) ■ (38) 

Within this class of models an important case arises for which the drift and diffusivity 
tensors are proportional to one another, 

km = TKlm ■ (39) 

In this case the microcurrent vanishes in equilibrium and the static probability distri- 
bution satisfies 

diKiUdi - Tdi(t>i^))Qoi^) = . (40) 

That is 

go(x) = iVe^*('') , (41) 
where N is a normalization constant. It is also the case that 

u^x) = / d^xr/€i^a„(/.(x)Qo(x) = . (42) 
Jn 

It follows that u(x) = u(x) . 



The formula for the effective diffusivity, eq(0), becomes 



kL = - X ^''^^''^'^'(^)^(^' x')n„(x')Qo(x') (43) 

which is symmetrical in I and m . Since G(x, x')(5o(x') is a positive operator we see 
that the tendency is for the random behaviour of the medium to reduce the effective 
diffusivity in a given direction relative to the microscopic value. This is consistent with 
the idea that the random fluctuations of the field 0(x) represent a trapping mechanism 
tending to hold the particle near the bottoms of potential wells. 
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5 External Gradient 



To investigate the computation of the effective drift we add a constant external gradient 
to the fluctuating gradient. We have 

ii;(x) = rKi„((9^0(x) . (44) 

The static probability distribution now satisfies 

diKim{dm - r((9m0(x) + gm))Qoi^) = . (45) 

When is small we can compute ui to lowest order in ■ If we set 

gg(x)=go(x) + gi(x) + --- , (46) 

then we easily find 

gi(x) = -y'dVG'(x,x')r/€;^a;Qo(x')^/m • (47) 
The mean velocity is, to 0{g), 

Ul=T (^Kim - J d^xd^x'TKindn(l){^)G{x,x')KrmdrQoi^')^ 9m ■ (48) 

From eq(^) we find 

= r ^Kzm - y c?'^X(i^x'ni(x)G'(x,x')Mm(x')(5o(x')^ 5fm ■ (49) 

This is also the result we obtained in a previous paper. It shows that the drift tensor 
is renormalized in the same way as the diffusivity tensor if we start from a situation 
where they are proportional at the microscopic level. Because this proportionality of 
the effective drift and diffusivity tensors holds sample by sample it also holds after 
averaging over the ensemble of random media. The precise nature of the statistics of 
the velocity field are not important for this result. There is no requirement for Gaussian 
statistics, for example. 

The result will also hold in the limit in which sample size becomes infinite, assuming 
the absence of very long range correlations. Consequently the result guarantees the 
Ward identity that was discussed in previous work |]T1|, [13], |T^. It remains remarkable 
that the direct proof of this identity, even in perturbation theory, is so unobvious that 
it has only been achieved at the two-loop level. It would be very interesting to obtain a 
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proof of the Ward identity directly in the infinite medium to all orders in perturbation 
theory. 

The results confirm those of Kravtsov, Lerner and Yudson 0, ^ in the case of isotropic 
diffusion, drift and (Gaussian) statistics. Our result shows that the only important 
point, as inferred by these authors, is the existence of the potential controlled sta- 
tionary distribution with a vanishing microcurrent. In fact by an appropriate linear 
transformation on the independent variables x, the diffusivity tensor can be rendered 
isotropic and simultaneously, the drift tensor rendered diagonal. Aside from the ques- 
tion of the statistics of the field 0(x) therefore, the more general problem considered 
here is mathematically identical to the isotropic case. 

6 Incompressible Flow 

Much work has been done on incompressible flow ^ ^. The velocity field satisfies 

diui{^) = . (50) 

The diffusion equation for the probability distribution can can be written 

|p(x.i) = (««a,-«.(x))a.P(x.«) . (51) 

It is then easy to see that the static distribution satisfies 

diQo{x) = . (52) 

This implies that Qq = 1/V . The role played by the mean velocity in this case is 
rather trivial. It is just the centre of mass motion of the incompressible fluid. For this 
reason, as noted in 0, it is never renormalized by the fluctuations in the velocity 
field. 

In an appropriate reference frame then, the mean velocity is zero. That is 

/ d^xu;(x) = . (53) 
Jn 

From eq(0) it follows that 

<m = ^^Im + rf^xrf^x'^iKx)G(x, ^')uU^')/V . (54) 

This is similar to eq(|13|) but with a change of sign for the correction to the bare 
diffusivity. This implies that the effect of the incompressible flow is to enhance the 
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diffusivity. In fact in practical cases the contribution of turbulent dispersion to the 
effective diffusivity dominates that of molecular diffusion by many orders of magnitude. 
In this case it is reasonable to consider a limit in which the molecular diffusivity goes 
to zero. This is not possible in the gradient flow case where the particle will simply 
become trapped at the first stationary point it encounters in the potential. 

7 Large Sample Limit 

For completeness we briefly discuss the large sample limit which we have assumed is 
innocuous in the sense that appropriate spatial averages are equivalent to ensemble 
averages. A simple example is the normalization of the static probability distribution 
in the case of gradient flow. We have 

Qo(x) = 7Ve^^(") . (55) 

Consider 

[ d^xe^-^W . (56) 
The ensemble average, in the case of Gaussian statistics, is 

(X) = / d^x(e-^(^)) = / (i^xe^'^(°) = T/e^'^(°) , (57) 

where A(x — x') = (0(x)0(x')) . We also have 

{X^) = [ (i^x(i^x'(e^<^W+^^(''')) = / ci^xci^x'e^'^(°)+^'^(^-^') . (58) 
Ja Jci 

The variance of X is therefore 

a\X) = I rf^xd^x'e^'^(°) L-^A(x-x') _ ^1 ^ y^r^HO) [ Lr^A(x) _ ^^1 _ .59) 

It follows that in the limit of large samples in which V ^ oo, both the mean and 
variance of X are OiV) . That is for large V we may use the result 

x = ye^^'^(°)(i + , (60) 

where ^ has zero mean and (7^(^) = 0{1/V), which suggests that ^ can be ignored in 
the large sample limit. This tells us that in this limit we can use the normalization 

Ar = e-i-'^W/y . (61) 
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A careful analysis of perturbation theory suggests that similar results hold for other 
appropriate quantities. In particular we expect for large samples, 



Mk) = (Mk)) + 0(l/v^) 



(62) 



Which implies that the infinite medium effective probability distribution is 



gik.X-(M(k)>t 



(63) 



This result was assumed in our discussion above. 

8 Spatially Varying DifFusivity 

So far we have applied the Derrida argument to models for which the diffusivity tensor 
is constant in space. Clearly the argument can be generalised to accommodate spatial 
fluctuations in the diffusivity. We therefore make the replacement 



and 1 is an element of the Bravais lattice. Of course we assume that the correlation 
length of the fluctuations in the diffusivity tensor is very much smaller than the size of 
the cell Q . For simplicity we will assume also that the local drift term vanishes since 
it will not be needed for the application we wish to make. The physical realizations 
of this problem are in the computation of the bulk versions of the diffusivity tensor, 
permeability tensor or permittivity tensor in the corresponding diffusion, fluid flow 
or electrodynamical problems 0] together with the bulk response to the appropriate 
external fields. 

If we adopt the diffusion model, the probability density function for a particle in the 
medium satisfies the equation 




(64) 



where 



fi:jj(x + 1) 



Kij (x) 



(65) 



d_ 

dt 



P(x,t) = diKij{x.)djP{^x,t) 



(66) 



It is immediately obvious that the static solution is 



Qo(x) = 1/V 



(67) 
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In this case the apphcation of the Derrida technique, along the hnes discussed above, 
leads us to consider the perturbation analysis of the equation 

i7(k)Q(x) = -/i(k)g(x) , (68) 

where 

Hik) = idi + ik.^Kij{x){dj+tkj) , (69) 

and ;u(k) is the eigenvalue that vanishes when k = . In the standard way we write 
/i(k) = /ii(k) + /i2(k) + ■ ■ ■ and Q(x) = Qo(x) + (5i(x) + ■ ■ ■ . However in the present 
case it is easily established that /ii(k) = and that 

Qi(x) =^A;i / rfVG(x,x')9;/s:i,(x')Qo , (70) 
where the Green's function ^(x, x') satisfies 

9ift:i,(x)6',G(x,x') = -5(x,x') . (71) 
A quick way to compute /i2(k) is to integrate eq(|68|) over Q to obtain 

/i(k) = -iki d^xKij{x){dj + ikj)Q{x.). (72) 
On picking out the 0{k'^) term from this equation, we find 

H20i-) = Kijkikj - iki J^d^xKij{-x.)djQi{-K) , (73) 



where the spatial average of the diffusivity tensor 

1 

V Jn 

From eq(^) we find 

1120^) = fiijhkj + kikj d^:>cd^x'Kii{x)diG{^,x')dl^Kjj^^ , (75) 
and after integration by parts 

/i2(k) = R,,hk, - ^ d^xci^x'[9,«,Kx)]G(x, ^')d'^K,U^') . (76) 
Hence we have obtained an expression for k^j. 

We now restrict ourselves to the isotropic case where = KQe^'^^^^dij. In the 

literature the case of lognormal statistics where 0(x) is a zero mean Gaussian field 
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with two point function A(|x — x'|) has been widely studied. It can also be related, 
by means of the appropriate Green's functions, to a diffusion problem with constant 
diffusivity kq and drift parameter Aq = tkq ■ The Green's function for this latter 
problem, G^™'^(x, x'), satisfies 

^M^^ - r9,0(x))G^™'^(x, x') = -5(x, x') , (77) 

whereas G(x, x') satisfies 

9i/€oe^'^W9iG'(x,x') = -5(x,x') . (78) 

It is easy to see now that the Green's function for the original diffusivity problem is 
related to that for the gradient problem by 

G(x, x') = e-"'^WG'3™'^(x, x'). (79) 

Substituting this into ([76|) we get 

/i2(k) = /tP - / c/^xrf^x'[9i/€oe"*(")]e-"^(")G^™'^(x,x')9'/€oe"'^("') , (80) 

where 

f^ = ^J rf^'x/toe'^^") . (81) 
If now we take the ensemble average k becomes 

« = i 1 (i^x/€o(e"^(")) = /€oe^^'^(°) , (82) 

and obtain the result 

= fi:oe^'^(°) - -r/t^Mi / d^xrf^x'([9i0(x)]G5™'^(x, x')9^e"'^("')) , (83) 
D V Jn 

We can compare this with the result of substituting = Kgdij and Mj(x) = rKo9j0(x) 
into eq{^. 

4m = f^oSim - r\l d^xd^x'[a;</.(x)]G'^™'^(x, x') [^^(/.(x')] Qo (x') , (84) 

where (5o(x) = Ne^'^^^^ . On taking the ensemble average and using the isotropy we 
obtain for the effective diffusivity in the gradient flow case 

^lra, = ^^-y'4^^^J^) /^d^xrfV[a„0(x)]G^™^(x,x')[9^0(x')]e^^(''') . (85) 
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The result is that 

= k^^^^/kq. (86) 

This means that the ratio of the bulk diffusivity to the (ensemble or spatial) average 
of the local diffusivity in the fluctuating diffusivity problem is equal to the ratio of the 
bulk diffusivity to the local diffusivity in the associated gradient drift problem. 
For the case of a lognormal diffusivity then we have 

^ e^^^(°)«:^,,,. (87) 

The consequences of this simple result are quite significant. In one dimension both 
problems are easily solved exactly and it is known that 

= Koe-"^-'^(^) (88) 

and 

^lrad = ^oe-'^^'^. (89) 

Hence our result is of course consistent with the well known exact results in one di- 
mension. 

In two dimensions however the random diffusivity case may be solved exactly via duality 



arguments |19, 20, 21], giving simply k° = 1. This means that n'trad ~ ^ 2^^^*^°^ giving 



^grad 

a new exact result for the isotropic gradient flow problem in two dimensions. 

In fact, these results in one and two dimensions are consistent with a renormalisation 



group calculation in D dimensions for the isotropic gradient flow case |13, 12| that gives 



^Irad = ^Oe--^^(°)/^. (90) 

Although the renormalization group calculation is based on certain plausible though 
unproved assumptions, it maintains the equality of the drift and diffusion renormaliza- 
tion factors and has been accurately verified in three dimension by numerical simulation 
over a considerable range in the value of this factor. It has been verified by analyt- 
ical calculation to two loop-order in perturbation theory. The numerical simulation 
suggests however that it remains accurate beyond this order. 

Eq(p6D allows us to relate these renormalization group results to the diffusivity problem. 
There has been for some time, a speculation in the literature, for example see ||18[, again 



verified to two-loop order in perturbation theory, that in the diffusivity problem 

«:« = «:oe(5-*)^'^W. (91) 
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Because eq(^) is an exact relation it shows that they stand or fall together. Recently 
2^ , p3| it has been shown that, in three dimensions, at third loop order eq(pl|) breaks 



down and corrections explicitly depending on the spatial structure of the correlation 
function appear. This would mean therefore that the formula ( pUD similarly fails at 
three loop order in three dimensions. So far it has not proved possible to analyse 
this problem to three loops directly. However a detailed analysis relating the two 
calculations should make further progress possible in the gradient drift problem. It 
is worth noting that even if a breakdown of the simple formulae is established, the 
success of eg (pOl) in describing the results of numerical simulations suggests that it is 
not necessarily a severe one. 

9 Conclusions 



We have developed the Derrida technique |T^, which was originally introduced for 



lattice models [jl], |T^, to the continuum case. It shows very conveniently the relation- 
ship between the effective drift and diffusivity tensors in general random media. In 
particular, for gradient flow models, it makes clear that when the microscopic drift 
and diffusivity tensors are proportional to one another the corresponding effective ten- 
sors show precisely the same proportionality [0. This guarantees the validity, under 
very general statistical regimes, of a Ward identity relating the effective propagator 
and vertex in wave vector space, derived up to the two-loop level in previous papers 
0, |15|, |16|. It ties in with the ideas behind the Einstein relation of the drift P, |^ 
and diffusivity by a temperature factor. Here the role of temperature is played by the 
inverse of the proportionality factor r . Of course this may have no direct connection 
with the real temperature of the sample. 

Applied to incompressible flow the method recovers immediately the non-renormalization 
of the mean flow by fluctuations in the velocity field. This is an obvious result on 
physical grounds and confirms previous work on the renormalization group approach 
to anomalous diffusion p|, ^. 

The method can also be applied to the problem of a random medium with a spatially 
fluctuating diffusivity. It reveals for the isotropic log-normal case the relationship 
of the bulk diffusivity with that for a related gradient flow model. This relation is 
confirmed in the case of one and two dimensional models and provides a constraint 
on the predictions of models in higher dimensions. The violation of simple results for 
the bulk diffusivity in one model being tied to possible failure of simple results in the 
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other. 

In principle the formulae derived from the Derrida method also allows the calculation 
of higher order terms that control the approach to asymptotic behaviour of the various 
moments of the probability distribution. However the exploitation of such results would 
require a proper elucidation of the commutation of the large sample and large time 
limits. This is an important mathematical problem that remains to be investigated. 
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